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ABSTRACT 

The accuracy of the numerical algorithm depends not only on the formal order of approximation 
but also on the distribution of grid points in the computational domain. Grid adaptation is a proce- 
dure which allows optimal grid redistribution as the solution progresses. It offers the prospect of 
accurate flow field simulations without the use of execessively time, computationally expensive, 
grid. 

Grid adaptive schemes are divided into two basic categories: differentiall and algebraic. The dif- 
ferential method is based on a variational approach where a function which contains a measure of 
grid smoothness, orthogonality and volume variation is minimized by using a variational principle. 
This approach provideds a solid mathematical basis for the adaptive method, but the Euler-La- 
grange equations must be solved in addition to the original governing equations. On the other hand, 
the algebraic method requires much less computational effort, but the grid may not be smooth. The 
algebraic techniques are based on devising an algorithm where the grid movement is governed by 
estimates of the local error in the numerical solution. This is achieved by requiring the points in the 
large error regions to attract other points and points in the low error region to repel other points. 

The development of a fast, efficient, and robust algebraic adaptive algorithm for structured flow 
simulation applications is presented. This development is accomplished in a three step process. The 
first step is to define an adaptive weighting mesh (distribution mesh) [1][2] on the basis of the equi- 
distribution law applied to the flow field solution. The second, and probably the most crucial step, 
is to redistribute grid points in the computational domain according to the aforementioned weighting 
mesh. The third and the last step is to reevalute the flow property by an appropriate search/interpolate 
scheme at the new grid locations. 

The adaptive weighting mesh provides the information on the desired concentration of points 
to the grid redistribution scheme. The evaluation of the weighting mesh is accomplished by utilizing 
the weight function representing the solution variation and the equidistribution law [3]. The selec- 
tion of the weight function plays a key role in grid adaptation [4] A new weight function utilizing 
a properly weighted boolean sum of various flowfield characteristics is defined. The redistribution 
scheme is developed utilizing Non-Uniform Rational B-Splines (NURBS) representation [5]. The 
application of NURBS representation results in a well distributed smooth grid by maintaining the 
fidelity of the geometry associated with boundary curves. Several algebraic methods are applied to 
smooth and/or nearly orthogonalize the grid lines. Elliptic solver is utilized to smooth the grid lines 
if there are grid crossing. 

Various computational examples of practical interest are presented to demonstrate the success 
of these methods. (1) Single wedge case. (2) Cylinder case. (3) Airfoil NACA 001 2 case. (4) Viscous 
missile case. (5) Slot cooled seeker window case. (6) Moving shock in converge inlet case. (7) 3-D 
single wedge case. (8) 3-D viscous cone case. (9) 3-D f 15 body and wing case. (10) 3-D helicopter 
propeller shock prediction case. In case (l),(2),(4),(5),(6),(7),(8),and(9), PARC2D andPARC3D [6] 
are used as the flow simulation code. In case (3), UBI [7] is used as the flow simulation code, and 
compare to Peter M. Goorjian’s solution [8], In case (10), MSUTC [9] is used as the flow simulation 
code. 
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Algorithm for Algebraic Adaptive Grid Generation System 







Figure 1. 

Use NURBS surface to generate smooth grid. 



Single Wedge Case 
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Figure 3. 

Single wedge grid adaptation process. 




Figure 4. 

Single wedge, closer Mach number contour look near shock region. 
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Figure 5. 

Use B -Spline curve to keep sharp corners. 
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Figure 6. 

Flow around cylinder initial/adaptive 
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Airfoil NACA0012 case n< 



Inviscid Pressure coefficients for NACA 0012 

1 

Mach-0.8, alphas 1 25, iir*3000, 120x100 Grid, i direction 
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Figure 9. 

Airfoil NACA0012 case Cp vs x/c plot, Mach = 0.8, a = 1.25. 



Viscous Missile Case 



Flow Solver Applied : PARC2D 








Figure 12. 

Slot cooled seeker window case, oblique shock above deflection comer. 
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Figure 13. 

Slot cooled seeker window case, velocity contour near interceptor. 
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Figure 14. 

Adaptive grid for the moving shock in the converge inlet. 
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Figure 15. 

3-D Single wedge, comparison of initial and adaptive grid & solution. 



3-D Viscous Cone Case 



Flow Solver Applied : FARC3D 
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Figure 16. 

3-D viscous cone case, bow shock in front of cone head. 


3-D F-15 Body and Wing 



Flow Solver Applied : PARC3D 
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Figure 17. 

3-D F-15 body and wing case, symmetric surface. 



Helicopter Propeller 
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